LTH-322 









> 



(N 

o 



X 



Anomalous dimension of non-singlet Wilson 
operators at 0{l/Nf) in deep inelastic scattering. 



J. A. Gracey, 
Department of Applied Mathematics and Theoretical Physics, 
C^ ■ University of Liverpool, 

~ P.O. Box 147, 

Liverpool, 

L69 3BX, 

United Kingdom. 



Abstract. We use the large A^ self consistency formalism to compute 
the 0{\/Nf) critical exponent corresponding to the renormalization of the 
flavour non-singlet twist two Wilson operators which arise in the operator 
product expansion of currents in deep inelastic processes. Expanding the d- 
dimensional expression in powers of e = (4 — d)/2 the coefficients of e agree 
with the known two loop structure of the corresponding renormalization 
group function and we deduce analytic expressions for all moments, n, at 
three and higher orders in perturbation theory in the MS scheme at 0{1/Nf). 



The quantum field theory describing the dynamics of strongly interacting 
particles is quantum chromodynamics, (QCD), which is an asymptotically 
free gauge theory, [|l|, 0]. In other words at increasingly higher energies the 
theory behaves more and more like a free theory and therefore the machinery 
of perturbative quantum field theory can be applied to gain an accurate 
understanding of physical processes. For example, the phenomenology of 
deep inelastic scattering can be explored and the one and two loop QCD 
corrections to parton model predictions can be computed which are found 
to be in good agreement with current experiments. (See, for example, |p.) 
With the increase in energy that will soon become available at LHC there 
would now appear to be a specific need to begin to examine subsequent three 
loop corrections to these processes, since the error bars on experimental 
results will be refined substantially. 

One of the central ingredients of the scattering formalism in gauge theo- 
ries is the operator product expansion, (OPE), of electromagnetic and other 
currents, |Q], which are composite operators of the fields of the standard 
model and QCD lagrangians. The decomposition of the Wilson OPE is in 
terms of an infinite sum of composite operators which are either flavour 
non-singlet or singlet, [H]. However, these operators require renormaliza- 
tion. In other words, radiative corrections generate a non-zero anomalous 
dimension, 7^8 (5)' ^°^ ^^^ ""^^ moment of, say, each non-singlet operator, 
where g is the coupling constant which will be defined later. Knowledge of 
7j^g((7) and the anomalous dimensions of the singlet operators are essential 
for determining, for example, corrections to sum rules. Although the Wilson 
coefficients are also required for this, where these are the coefficients which 
appear with each operator of the OPE, they depend on the specific process 
involved, unlike the anomalous dimensions, and therefore require separate 
treatment, |p. Like the functions of the renormalization group equation, 

(RGE), 7ns (5) ^^^ calculated order by order in perturbation theory. Specif- 
ically one renormalizes the operator as a zero momentum insertion in some 
Green's function. As a result the one loop structure of non-singlet opera- 
tors, 7j«jq(5'), was given in 0] whilst the two loop analysis was carried out 
in [5-7]. (Two loop results for the singlet functions have been deduced in 
[4,8-11].) So it is the perturbative structure of these functions which must 
be determined at three loops in order to refine the theoretical understanding 
of physical processes. 

Such calculations, however, are formidable in their computational com- 
plexity due to the huge number of Feynman diagrams to be analysed. In- 



deed the two loop results were performed using computer algebra packages. 
Therefore it is also important to have an independent way of calculating in- 
formation on the perturbation series which agrees with known calculations 
and, equally, provides additional non-trivial information at much higher or- 
ders. In this letter we introduce such a technique for computing the anoma- 
lous dimensions of physical operators involved in the OPE of deep inelastic 
scattering, though for simplicity we will only concentrate on the non-singlet 
case to illustrate the power and beauty of the method. (As the Wilson co- 
efficients are process dependent they are beyond the scope of this letter.) 
The formalism is based on the large A^ self consistency programme which 



was first introduced in |12, 13 1 for the 0{N) bosonic a model. To deduce 
information on the perturbative structure of any function of the RGE, one 
computes the relevant critical exponent at the non-trivial zero, gc, of the d- 
dimensional /3-function of QCD. There the theory is finite and the propaga- 
tors of the fields and Green's functions have a simple power law or conformal 
structure. The associated critical exponent is, by universality, a function of 
the spacetime dimension, d, and any internal parameters. For the present 
letter, these will be A^, the number of quark flavours, and Ac the number of 
colours. The benefit of proceeding in this fashion is that the RGE simplifies 
because of P{gc) = 0, Oc ¥" 0; to the extent that the d-dimensional critical 
exponent is simply related to the corresponding RGE function evaluated at 
the critical coupling. (See, for example, [^4[.) Therefore, if one computes 
in some approximation, which will be large A^ here, then knowledge of the 
critical exponent in d-dimensions means calculating order by order in 1/Nf 
one can deduce the coefficients of the perturbative function away from crit- 
icality. At low orders of g the coefficients deduced from the exponent will 
agree with those calculated in the mass independent MS scheme which uses 
dimensional regularization. Clearly this will be a powerful way of proceeding 
since one will always be able to obtain information to all orders in g at each 
level in 1/Nf. Indeed it is one aim of this letter to draw attention to this 
alternative method of calculating the perturbation series of useful physical 
quantities. 

We recall for the interested reader that the earlier applications of the 
large A^ exponent programme to four dimensional gauge theories included 
the evaluation of the electron anomalous dimension at 0{1/N?), [^], the 
QED /3- function at 0{1/Nf), [|l^], and the electron mass anomalous dimen- 
sion at 0(1 /N?), ^^. As the latter involved the computation of the expo- 
nent of the mass operator ?/J^ which is also composite, we have used [^] 



as a foundation for the present work. More recently, the quark, gluon and 



ghost anomalous dimensions have been deduced at 0{1/Nf) in |15] which 
are in agreement with the previous three loop results of [1,19-22]. 

To proceed with the critical point approach we recall the form of the 
QCD lagrangian we use is 

L = - -^^^ + if^q'^ + AlfTfj-f^'q'-^ - \ f^'^daAlA''^ A"^ 

rabc fade Ab \c Afid/tue / q Afj,a\2 

- d^c^d^.c" + ^'"'^^,c^c^A^''' (1) 

where g*^ is the quark field, I < i < Nf, 1 < I, J < N^, A^^^ is the gluon field 
with I < a < N^ - I, F^^ = d^Al - d^A^, & and c" are the ghost fields, 
^ is the covariant gauge parameter, Tfj are the generators of SU{Nc) which 
has structure constants f"'^'^ and e is the coupling constant. The location of 
Qc in the neighbourhood of which we will base our analysis is given by f3{gc) 
= 0, where to one loop in d-dimensions, M, 



(3{g) = {d-4)g 



lTiR)Nf - ^C2iG) 



9^ + 0{g^) (2) 



where the dimensionless coupling is g = (e/27r)^, in the notation of [^], and 

Tr(r'*r'') = T{R)5°'^ , T^T"" = C2{R)I , f''^'^ f"^ = 02(0)5"^ (3) 

are the Casimirs of a general classical Lie group. For SU{Nc), T{R) = i, 
C2{R) = (iV2 - l)/2iVc and CaCG) = N^. Therefore from (2) 

where d = 4 — 2e but it turns out that to deduce the coefficients for 7j;^g (g) 

from the exponent r/^"' = j^Mgc) we calculate at 0{1/Nf) here, the 0{1/N?) 
corrections of (4) are not required. 

The Wilson operators whose critical exponent r/'"' we will compute are 
the twist 2 flavour non-singlet operators, |Q], 

^NsT'^ = 5^""W7^'^'''---^''"7Tj(1 ±7^)9-^ -trace terms (5) 

where D^ = d'^ + iT"'A^^"' and S denotes the symmetrization of the Lorentz 
indices. At leading order in large Nj it will turn out that the 7^ term 



will be completely passive and the same result obtained for either sign. So 
we will suppress the index it for the moment. To determine the anoma- 
lous dimensions one first of all removes the symmetrization and trace terms 
by introducing a null vector A'^ with A^ = and multiplies Oj^g'^" by 
A^^ . . . A^^ to simplify the number of insertions to be made in the Green's 

function, (^Ons^! 0- Next as the anomalous dimension 7^3(5) is gauge 
independent one chooses to compute in the Feynman gauge to reduce the 
algebra involved in the calculation. As a consequence of both these steps 
the insertion of the operator (5) requires the development of new Feynman 
rules. These have been derived in Q for the one loop case and in ^ for the 
2-loop calculations. We will use the former rules for the large Nf computa- 
tion since at leading order the only graphs which contribute are those of fig. 
1, where the dotted line corresponds to quark fields. The circle with cross 
denotes the zero momentum insertion of Ons- 0-^ the singlet calculation 
there would be additionally two 2-loop graphs with a closed quark loop, 
[^, contributing at this order in l/Nj, but these are zero in the non-singlet 
case.) 

In the perturbative calculation one computes each graph of fig. 1 with 
the conventional quark and Feynman gauge gluon propagators, ^//c^ and 
'Hfj-u/k'^ respectively, in dimensional regularization and absorbs the simple 
poles in e minimally into the appropriate renormalization constant. In the 
critical point approach one uses an alternative strategy. Near g^, (4), the 
propagators of the fields obey asymptotic scaling and therefore satisfy a 
simple power law structure consistent with Lorentz and conformal symme- 



try. The general forms have been used extensively before, |15, 18], and in 



momentum space we recall that they are, in the Feynman gauge, 

where A and B are momentum independent amplitudes, ji = ^d and the 
critical exponents a and f3 are defined as 

a = ix-l + \r] , fi = l-n-x (7) 

Here the canonical dimension of each field is deduced from a dimensional 
analysis of the kinetic terms and interactions of (1) and the anomalous terms 
r/ and x ^-re respectively the quark wave function renormalization exponent 
and the qqg vertex anomalous dimension. (At leading order in 1/A^ there are 
no graphs involving ghosts, [H].) The presence of the non-zero anomalous 



dimensions means that when considering subsequent orders in 1/A^ only 
those graphs where the quark and gluon propagators are not dressed are 
considered. From |jl^, O] we recah that in the Feynman gauge 

"" {2^,-l)i^i-2)T{R) ^^^ 

where r? = E^=i Vi/N} and 7?° = {2fi - 1) (/x - 2)r{2fi)/[iT^ {fi)T{fi + l)r(2 - 
/i)]. Also, the combination z = A^B is required and [|l^, |l^ 

' 2{f,-2){2f,-l)T{R) ^' 

We now recall the method to compute the anomalous dimensions of 
composite operators with the critical propagators (6). First, one substitutes 
the lines of fig. 1 with (6) and calculates the integral. However, one would 
discover that the integral is divergent and thus requires regularization. This 
is achieved by shifting the gluon exponent by an infinitesimal quantity 6, 



f3 —^ P — 5, |l^, |T^, 18]. (In previous work, we used the symbol A for 
the regularizing parameter but to avoid confusion with A^ we use 5 here.) 
When one computes with non-zero 6 now, each graph of fig. 1 will have 
the following formal structure at 0{1/Nf) after expanding in powers of 6, 



@0' 

^ + Q + Rlnp^ + 0{5) (10) 



where P, Q and R depend on ^ and Nc, and p is the external momentum 

flowing through the quark fields. The simple pole is absorbed minimally by a 

conventional (critical point) renormalization, |24|. This leaves a term involv- 



ing R which would violate scaling symmetry in the limit as one approaches 
criticality, p^ — > oo. To avoid this one notes that by general arguments 
the (5-finite Green's function must resum to the structure {p'^)'^o (9c)/^^ p^ , 
where Jq (gc) is related to the exponent we are interested in, rji , through. 



m 



ryW = r? + 7&^(5c) (11) 

and is what is obtained from the insertion of (5) in the graphs of fig. 1. 
Therefore, by isolating and then summing the contributions from the Inp^ 
terms of each graph of fig. 1, we can deduce r/}" via (8) and (11). 



Using (6) and the Feynnian rules of Q, we find that the first graph of 
fig. 1 contributes 



2M/x-l)=^C2(i?K 



in -2){2fi-l){n + n-l){n + n- 2)T{R) 



(12) 



.(") 



to 70 (5c)) whilst the second (together with its mirror image) gives 



ifi + l 



{^x-2){2^i-l)T{R) ^^^ 
Therefore, from (8) and (11) we have 

2C2{R){^l - Ifrii \ (n-l)(2/i + n-2) 



(13) 



(n) 



(2/i-l)(/.-2)T(i?) 



+ 



2n 



(/^ 



n^ 



{fi + n- l)(/i + n 
2fi 



f^_^{^^ + l-2) {is-lf 



(14) 



which is the main result of this letteiQ. 

There are several checks we can make on (14). First, in the case n = 1, 
the anomalous dimension 7^3(5) is zero, [Q], since then O^g corresponds to 
a conserved current. It is easy to verify from (14) that 77]^ = in agreement 
with this general result. Second, we can expand r/J"' in powers of e in d = 
4 — 2e dimensions and compare with the 2-loop result of [4-7] ie 



(n),± 



7ns (5) 



C2{R)g 



2Siin) 



1 



2 n(n + l) 



C2{R)T{R)Nf 



3 

+ 4 + 



;iln2 + 5n-3)^ 
rfl{n + 1)^ 



652(n)-10Si(n) 



+ b^{n) 



(15) 



where Si{n) = Y^^=i 1/^' a-^id the coefficients b^{n) have been given in [5-7] 
but are 0{1/N?) and their explicit form is not needed for the present point. 
Substituting the critical coupling gc in (15) from (4) and comparing with 

^Although the factor of (/i — 2) in the denominator of (14) suggests that r/j behaves 
as 1/e as e ^ when /i = 2 — e, in fact r]i = 0{e^) which means that 77^ and J-i^l^igc) 
are both 0(e). 



the e expansion of (14) it is easy to verify agreement with the two 0{1/Nf) 
coefficients of (15). This is a highly non-trivial check and establishes the 
correctness of (14) since we considered only the leading order 1/Nf graphs 
which are one loop. In other words the conformal ansatze and formalism 
correctly reproduce the large A^ bubble sum contributions at higher order. 
More importantly having verified its correctness we can now deduce subse- 
quent terms in the perturbation series albeit at leading order in 1/Nf. For 
example, if we set 



^^^>{g) = aiC2{R)g + Y.a,C2{R)[T{R)NfX-^g' 



(16) 



for the leading order 1/A^/ part of the function, where ai and 02 can be read 
off from (15), then (14) implies that 



03 



04 



ag 



2 10 

-S'3(n) - i^S2{n) 



2 17 

27^^^")+ 72 



[12n'^ + 2n^ - I2n'^ - 2n + 3] 



27n3(n + l)2 



Si{n) 



10 



^S,(n) 



[4n« 



27 
+ 

+ 



^[355(n 
C(4) 



81 



S.(n) 



81 



-S2(n) 

1 
9 



81 



Si{n) 



C(3) + 



27n(n + 1) 
12n^ - 15n^ + lOn^ + 14^^ 



131 
1296 

n-3] 



81n4(n + l)4 
554(71.) - 53(72) - S2{n) 



Si{n)] 



+ 



+ 



+ 



54 

03) 
486 

[87Z^ 



45i(7i) 



2452 (ri) 



n{n + 1 
405i(n)+3 + 



+ 



323 



7776 
4[lln2 + 572-3] 



Ti?{n + 1)^ 
Sn^ - 28n5 - 571"^ + 2477,^ + 1377,^ - An - 



2437i5(n + l) 



(17) 



(18) 



(19) 



are the subsequent new coefficients, where C,{q) is the Riemann zeta func- 
tion. Of course, 03 will be an important check for the explicit 3-loop MS 
calculation of 7^0(9) which is currently being calculated for various mo- 
ments, [^]. To aid comparison of the 0{1/Nf) part of the full three loop 
result of [p5|], we have evaluated 03 for even ?i, 2 < ?i < 22, and collected the 



results in table 1. In p5|], the coefficient of the n = 8 case was given and it is 
very satisfying to record that there is exact agreement with the correspond- 
ing entry of our table, which reinforces our confidence in the correctness of 
the exponent (14), (after allowing for the different convention in defining 
the coupling constant in this paper and an overall factor of 2 between the 
definition of the renormalization group functions) . The remaining values of 
table 1 will occur in the explicit evaluation of 7)00(5) in perturbation theory 
for other n. 

Finally, with (14) the x-behaviour of the 0{\/Nf) part of the non-singlet 
splitting function, P]>,j-g(g,a;), of the Altarelli-Parisi equations can now be 
examined, where x is the Bjorken scaling variable. This will be useful in the 
extension of the earlier 2-loop work of [0, 26 1. It is defined as the Mellin 
transform of the anomalous dimension, 

fdxx^-^P^^{g,x) = - hNs(5) (20) 



where we use the conventions of 127 



We conclude by noting that we have given an alternative way of com- 
puting anomalous dimensions of physical operators in QCD, which agrees 
with low order results, and also provided the means to deduce the higher 
order structure. It ought now to be possible to develop this large A^ formal- 
ism in two directions. First, the 0{1/N?) corrections can, in principle, be 
deduced. However, one would first need to know r/2, the quark anomalous 
dimension, which has yet to be computed. Secondly, it may be possible 
to determine the singlet operator anomalous dimensions. Although these 
mix under renormalization that is not a significant problem since matrices 
of exponents have been deduced in other contexts within the critical point 



programme, |24]. 
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n 


as 


2 


28 
243 


4 


384277 
1944000 


6 


80347571 
333396000 


8 


38920977797 
144027072000 


10 


27995901056887 
95850016416000 


12 


65155853387858071 
210582486065952000 


14 


68167166257767019 
210582486065952000 


16 


5559466349834573157251 


16553468064672354816000 


18 


19664013779117250232266617 


567701 18727793840841472000 


20 


6730392290450520870012467 


18923372909264613613824000 


22 


16759806821032136669044226177 


46048135637404510767879321600 



Table 1: Values of a^ for various n. 
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Figure Captions. 

Fig. 1. Leading order graphs for r?}"" . 
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